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Abstract 

A sizable rate of events where two pairs of 6-quarks are produced contem- 
porarily is foreseen at the CERN LHC, as a consequence of the large parton lu- 
minosity. At very high energies both single and the double parton scatterings 
contribute to the process, the latter mechanisms, although power suppressed, 
giving the dominant contribution to the integrated cross section. 
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I. INTRODUCTION 



One of the main topics at the LHC is the production of 6-quarks, both to search for 
CP violation, looking at b decays, and to test QCD, by studying the production mechanism 
Bottom quarks are also a large source of background to several processes of interest, as 
in various promising channels for Higgs detection [0]. The production mechanism of heavy 
quarks in hadronic collisions pose, on the other hand, non trivial problems already at smaller 
energies, where also the simplest observable quantity, the integrated inclusive cross section, 
is not reproduced trivially. 

The inclusive cross section of 6-quarks production has been evaluated in pQCD at the 
next to leading order in 0:5 0]. Unfortunately comparisons with the recent experimental data 
of the DO Collaboration H at TEVATRON have shown that the NLO pQCD calculations 
underestimate the cross section by a factor ~ 2, 3, showing that NNLO corrections, whose 
explicit evaluation is still an open question, give a large contribution to the cross section. 
A complementary approach to heavy quarks production, which keeps explicitly into account 
that transverse momenta and virtualities of the interacting partons become increasingly 
important in the kinematical regime of s ^ ml ~ s ^ A^, and includes terms at every order 
in as in the calculation of the cross section, is the fc^-factorization, where the interaction 
is factorized into un- integrated structure functions and off shell matrix elements |^ [Q. 
From the phenomenological point of view A;rfactorization is not inconsistent with HERA and 
TEVATRON data, allowing one to reproduce both the value of the integrated inclusive cross 
section and various differential distributions, including the correlation in the azimuthal angle 
between the produced b quarks, where different approaches are less successfully compared 
with experiment (see and references therein). 

Interestingly, although the value of the integrated inclusive cross section cannot be ob- 
tained trivially, one may find several cases where the overall effect of higher order corrections 
amounts to a simple rescaling of the lowest order parton model result. There are in fact 
several distributions, derived either using the /ct-factorization approach or by working out 
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the cross section at the NLO pQCD, which are rather similar (apart form the normahzation 
factor) to those obtained with a simplest lowest order calculation [^]. Hence, in a few cases, 
the whole effect of higher order corrections is (approximately) reduced to a single numerical 
value, the so-called K factor: 

= (1) 

where a(bb) is the inclusive cross section for bb production and aioipb) the result of the 
lowest order calculation in pQCD. 

When looking at extrapolations of the cross section at high energies, one finds that the 
result is affected by several uncertainties, as the knowledge of the parton structure functions 
at very small x and the values of the heavy quark mass and of the running coupling constant. 
Although the expected inclusive cross section of h production is hence still pretty uncertain 
at LHC energy, all estimates point in the direction of rather large values, as a consequence 
of the high parton luminosity The fairly large fiux of partons make it also plausible 
to expect a sizable rate of events, where two or more hh pairs are produced contemporarily 
by different partonic collisions in a given pp interaction []TU[. Although at present stage all 



quantitative predictions for this much more structured interaction process are unavoidably 
pretty uncertain, the large cross sections foreseen at the LHC is, in our opinion, a strong 
motivation to make an attempt of giving a few quantitative indications on the production 
rate of multiple hh pairs trough multiparton interactions at the LHC, comparing with the 
rates to be expected by the more conventional single parton scattering mechanism. 

Since the details of the elementary production of heavy quarks are still a matter of 
debate, we limit our considerations, for the production of multiple hh pairs, to the simplest 
cases, where the whole effect of higher order corrections is taken into account by the overall 
normalization factor. Given the lack of information on higher order corrections in the 2^4 
processes, we make moreover the assumption that the K factors of the gg —>■ hbhb and 
of the gg hh processes are equal. Hence we work out the gg —* hh process in the kt- 
factorization approach, fixing the input parameters by comparing with the TEVATRON 



data, and extrapolate the cross section at LHC energies, identifying a few distributions 
where the effect of higher order corrections reduceds to a simple rescaling of the lowest 
order result. The value of the A'-f actor derived in this way is then used to renormalize 
the double {gg —>■ bb)^ and the single gg —>■ bbbb parton scattering cross sections, which we 
evaluate by working out all Feynman diagrams at order a^. 

II. BB CROSS SECTION AT TEVATRON AND LHC AND K-FACTOR 

In the fct-factorization approach the bb production cross section is expressed as [0 

a{pp-^bb) = [ ^^-^^^-^dxidx2f{xi,qti,fx)f{x2,qt2,l^)cr{xi,qti;X2,qt2;i^) (2) 
J n n 

where f{x,qt,^) is the unintegrated structure function, representing the probability to find 
a parton with momentum fraction x, transverse momentum qt at the factorization scale /j,, 
while a is the off-shell partonic cross section of the process g*g* —>■ QQ. 

To work out the inclusive cross section we use two different prescriptions for constructing 
the fcf-distributions from the usual integrated parton densities. The first prescription is based 



on the conventional DGLAP evolutions equations ||rT|, with virtual corrections re-summed 
in the survival probability factor Ta{k'^,fi'^) [jl2|. Hence the un-integrated structure function 
for the parton a reads 



fa{x,kt,^ ) = Ta{kt,^ ) 



(3) 
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where Paa'{z) is the splitting function, a'{x,k^) the integrated structure function and 6 a 
cutoff parameter introduced to give sense to the integral. Although not written explicitly 
also Taik'j^n'^) depends on 5, in such a way the fa^Xjk'f,^'^) is a smooth function when 5 
becomes small. 

As for the second prescription we follow ref. [|T3|, where the un-integrated structure func- 



tions are obtained from the leading order BFKL equation and are expressed as the convo- 
lution of the usual colhnear gluon densities G{x, iJ?) with the universal function Q{x, k^, fj,"^) 
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H^,kl^^^) = C dig{i,kl^i^)G{^,^x^) (4) 

J X X 

The weight factors Ql^jkfjfi"^) have a known analytic expression, in double-logarithmic 
approximation, in terms of Bessel functions and depend on the quantity = SOs/vr, which 
in the BFKL formalism is a fixed parameter, related to the pomeron intercept a(0) = 1 + A, 



where A in leading log approximation is A = 4a5log2. Following [|I4[ we take the value 
A = 0.35. 

To generate the un-integrated structure functions we have used the parton distributions 
set GRV94 with factorization scale = s. Hence in evaluating the cross sections, in the 



fct-factorization approach, we have set the renormalization scale equal to the gluon virtuality. 
To obtain the cross section at the lowest order in pQCD we have used the MRS99 parton 
distributions , with factorization and renormalization scale equal to the transverse mass 
of the 6-quark. Comparing the total cross section values in the two approaches, we have 
obtained for the i^-factor the value K ~ 5.5. 

In Fig.|l] we plot the integrated cross section of bb production at TEVATRON(^ = 
l.STeV) and at LHC(y^ = 14TeV), as a function of the minimum value of the transverse 
momentum of the 6-quark. The dotted curves represent the cross section derived 
using the unintegrated gluon structure function, according with the BFKL prescription of 
Eq(4), whereas the dashed lines have been obtained by using the prescription in Eq(3). The 
continuous lines represent the result of the lowest order calculation multiplied by the K 
factor. At TEVATRON energy the 6-quark distributions are within the rapidity interval 
|?/| < 1 and are compared with the DO experimental data [Q. The same distributions, 
extrapolated at LHC energy, are then plotted as a function of within the pseudorapidity 
interval \ri\ < 0.9, corresponding to the acceptance of the ALICE detector. 

In Fig.|^ we show the rapidity(?/) and pseudorapidity(r7) distributions normalized to one 
and within |?7| < 0.9. Here the continuous histograms refer to the result of the lowest 
order calculation, rescaled by the K factor, whereas the dashed histograms represent the 
distributions evaluated with the fcrfactorization approach. As one may see also in this case 



the whole effect of higher orders reduces to a simple rescaling. 



III. BBBB CROSS SECTION 



The leading order QCD process to produce two pairs of heavy quarks is given by the 



single parton scattering term at the fourth order in the coupling constant ||T^. A com 



peting mechanism at the LHC energy is the double parton scattering [|18[. We compare 
the two mechanisms in proton-proton collisions in the kinematical range of the ALICE and 
of the LHCb detectors, namely at center-of-mass energies of 5.5 and 14 TeV, within the 
pseudorapidity regions |?7| < 0.9 and 1.8 < rj < 4.9, down to very low transverse momenta. 

The single scattering pQCD sub-processes at the lowest order in a.^, in pp —>■ bbbb, are 
the quarks initiated process, qq bbbb, whose amplitude is given by the sum of 14 Feynman 
diagrams for each flavor in the initial state, and gluon fusion, gg bbbb, represented by 
76 diagrams altogether, the latter amplitude giving the dominant contribution to the cross 
section at small x. To evaluate the cross section we have generated the matrix elements 



of the partonic amplitudes with MadGraph and HELAS |^ and we have used the 



MRS99 parton distributions |T6|, with the factorization scale equal to the renormalization 
scale fip = fJ'B., which we have kept fixed at the value of the transverse mass of the produced 
b quark. For the mass of the bottom quark we have used the value rrib = 4.6 GeV. The 
multi-dimensional integrations have been performed by VEGAS [0] and the resulting cross 
section has been finally multiplied by the K factor obtained as described in the previous 
section. 

The evaluation of the double parton scattering contribution to the cross section is con- 
siderably more uncertain because of the unknown non perturbative input to the process, 
given by the two-body parton distribution functions r{xi,X2, (3) |T^, where Xi^2 are the 
fractional momenta of the two partons belonging to the same hadron and (3 their distance 
in transverse space. Although not explicitly written, the distributions depend also on the 
scale factors characterizing each elementary interaction and on the different kinds of partons 
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involved. Given the large parton population at low x, to proceed further we make the usual 
simplifying assumption of neglecting correlations in fractional momenta and we factorize the 
two-body parton distribution as T{xi,X2, (3) = G{xi)G{x2)F{P), where G{x) are the usual 
one-body parton distributions and F{l3) is a function normalized to 1 and representing the 
parton pair density in transverse space. With these assumptions the cross section acquires 
the simplified form ^2 



aDihh- hh) = W &^a,{hh)a,{hh) (5) 

where the indices i,j label the different cases where each hh pair is originated either by a 
qq annihilation, discriminating the cases of sea and valence, or by two gluons and (Ji{bb) 
represents the inclusive cross sections for bb production in a hadronic collision, with the 
index i labelling a definite parton process. The weight factors O*-' have dimension an inverse 
cross section and result from integrating the product of the two-body parton distributions 
in transverse space, while the factor 1/2 is a consequence of the symmetry of the expression 
for exchanging i and j. The dependence of B*-' on the indices i,j accounts for the possibility, 
for different pairs of partons in the hadron structure, to be characterized by different values 



of their relative average transverse distance [|T^]. Notice that by measuring the double 



parton collisions one has access to a new information on the hadron structure, summarized 
in these weight factors, which cannot be obtained in hard processes with a single parton 
interaction only. 

The experimental information on the double parton scatterings is due to the four-jet pro- 
duction measurement in pp collisions at ^/s = 63 GeV, performed by the AFS Collaboration 
p3| , and to the study of final states with three minijets and one photon in pp collisions at 
y/s = 1800 GeV, due to CDF In both cases the cross section was expressed as 



mas{A)crs{B) 
(Jd = — (6) 

where m = 1 if the two parton processes A and B are identical, while m = 2 if they are 
different and as is the single scattering inclusive cross section. The overall output of the 



experiment hence reduces to the value of a single parameter, the scale factor (Jeff, whose 
value is aeff = 5 mb for AFS, while a^ff = 14.5 mb for CDF. The two experimental results 
are not inconsistent, given the different content of partons in the two cases, mainly valence 
quarks in the former case and mostly gluons and sea quarks in the latter; the experimental 
indication hence pointing in the direction of a sizable dependence of the factors G*-' in (j^) on 
the different elementary processes. Interestingly, the measurement of de// for different final 
states, as a function of the cm. energy and cuts applied, allows one obtaining the values 
of the scale factors G*-', letting in this way access to the three dimensional structure of the 



proton |18 



The dominant contribution to bbbb production at the LHC is gluon fusion, so, for the 
present purposes, the sum in (^) may be approximated well by a single term, where the 
scale factor could not be strongly different with respect to the case of the CDF experiment. 
Hence, to evaluate the double scattering cross section, we have used the simplest expression 

where for aeff we have taken the value reported by CDF. Notice that since an is propor- 
tional to (t|, the effect of higher order corrections is enhanced on au: 



as = Kal"" 

an = a^rF (8) 
where ag'^ refers to the lowest order expressions of the cross section. 

IV. RESULTS 

In Fig.^ we plot the expected rise of the total bbbb production cross section as a function 
of the cm. energy. The continuous curves refer to the double parton scattering contribution, 
while the dotted curves to single scattering. In each case the lower curve refers to the value 
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K = 2.5, while the higher curve to = 5.5, which are the typical estimate of the NLO-QCD 
and the result of our calculation within the fc^-factorization approach. Notice that at the 
LHC the double parton scattering gives a contribution to the integrated cross section about 
ten times larger than the single scattering. 

To see how the cross section depends on the transverse momenta we have plotted in Fig.^ 
the two contributions to the integrated cross section, as a function of p™*", the minimum 
value of the transverse momenta of the b quarks (which we require to be all inside the 
pseudorapidity interval \ri\ < 0.9), for the center-of-mass energy values of 14 and 5.5 TeV. 
The continuous histograms refer to the double parton scattering contribution, while the 
dotted histograms to single scattering. The double parton cross section decreases faster 
with p™" than the single parton cross section, the two contributions being of the same order 
at p™" = 8 — 10 GeV. Pseudorapidity, and rapidity distributions at 14 and 5.5 TeV are 
plotted in Fig.^ ( always requiring for the two 6-quarks \ri\ < 0.9), where continuous and 
dashed histograms have the same meaning as in the previous cases. 

To see how the results depend on rapidity, we have plotted in Figj^ the same rapidity (y) 
and pseudorapidity (77) distributions at -^i = 14 TeV, requiring both 6-quarks to be in the 
pseudorapidity interval 1.8 < t] < 4.9, which corresponds to the acceptance of the LHCb 
experiment. In the same figure we also compare the two contributions of single and double 
parton scattering, integrated within the rapidity acceptance of the LHCb, as a function of 

Pt 

The overall indication which one obtains from the present study is that double parton 
scatterings dominate the bbbb integrated cross section by a large factor, both in the central 
rapidity region and at the larger rapidity values of the LHCb experiment. In both cases 
the contribution of the single parton scattering term becomes important only after applying 
cuts to the transverse momenta of the order of 8-10 GeV. The large values expected for the 
bbbb cross section, which at 14 TeV are of the order of one /ib, inside the ALICE and the 
LHCb detectors, and the localization of the double scattering contribution at relatively low 
Pt values, suggest that heavy quark pairs production at the LHC might represent an efficient 
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tool for studying the gluon initiated double parton scattering process. 
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FIGURES 
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FIG. 1. pp ^ bb production cross section as a function of p™™ at y/s = 1.8 TeV, with the 6-quark 
within the rapidity range \yb\ < 1, experimental data from ref. and at y/s = 14 TeV with the 6-quark 
within the pseudo-rapidity range \r]\ < 0.9 . 
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FIG. 2. Normahzed rapidity (y) and pseudorapidity (rj) distributions for bb production at ALICE. 
With the kt factorization approach (dashed histograms) and at the lowest order in pQCD multiplied by the 
if -factor (continuous histograms). 
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FIG. 4. bbbb production cross section at ^/s = 14 TeV and at -^i = 5.5 TeV as a function of pj'"" with 
all the four 6-quarks in the pseudo-rapidity interval \r]\ < 0.9 . 
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FIG. 5. bbbb production with the two equal sign 6-quarks in the pseudo-rapidity interval \r]b\ < 0.9. 
//-distributions and y(,-distributions at ^/s = 14 TeV and at ^/s = 5.5 TeV. The continuous histograms refer 
to the contribution of double parton scatterings while the dashed histograms to the single parton scatterings. 
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FIG. 6. bbbb production with the two equal sign 6-quarks in the pseudo-rapidity interval 1.8 < 77 < 4.9 at 
^/s = 14 TeV. Production cross section as a function of p™'", r] and y. The continuous lines and histograms 
refer to the contribution of double parton scatterings while the dashed lines and histograms to the single 
parton scatterings. 
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